A dilute atomic Fermi system with a large positive scattering length 
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We show that a dilute atomic Fermi system at sufficiently low temperatures, can display fermionic 
superfluidity, even in the case of a repulsive atom-atom interaction, when the scattering length is 
positive. The attraction leading to the formation of Cooper pairs is provided by the exchange of 
Bogoliubov phonons if a fraction of the atoms form a BEC of weakly bound molecules. 

PACS numbers: PACS numbers: 03.75.Ss, 03.75.Mn, 03.75.Hh 
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Ever since the pioneering experiment of DeMarco and 
Jin there has been an ongoing search for the super- 
fluidity in a dilute fermionic atomic system . A recent 
experimental result from Jin's group Q seem to point 
to the possible creation of an even more intriguing sys- 
tem, a dilute atomic Bose system, subsequently shown 
to be in a BEC state H, coexisting with a dilute atomic 
fermionic fluid, perhaps even a superfluid itself as we 
will show here, both supcrfluids composed from one sin- 
gle atomic species. Taking advantage of a Feshbach res- 
onance between given hyperfine states all these authors 
0,0] have manipulated the scattering length and made it 
positive. Starting with an almost equal mixture of these 
hyperfine spin sub-states and a relatively small and neg- 
ative scattering length, by ramping the magnetic field at 
various rates, these authors have created a mixture of 
weakly bound fermionic atoms and unbound fermions, 
all of them interacting with positive scattering lengths, 
as we shall show here. The weakly bound state formed 
is bosonic in character, with a binding energy typically 
larger than the temperature of the atomic cloud |4] . This 
binding energy is approximately given by ti 2 /ma 2 , when 
the scattering length is large. We shall often refer to these 
weakly bound molecule as dimers. The fraction of such 
molecules depends on the rate at which the magnetic field 
is ramped across the Feshbach resonance, ranging from 
essentially zero in the sudden limit and up to about 80% 
in the adiabatic limit. At T = and in the adiabatic 
limit one would expect a full conversion of all fermions 
into such weakly bound molecules, the so called BCS 
— * BEC crossover [S IS HI 0] • Unless the magnetic 
field is very close to the Feshbach resonance, the size of 
these molecules is smaller than the average inter particle 
separation. If the weakly bound molecules form a BEC 
state, then the induced fermion-fermion interaction, due 
to the exchange of Bogoliubov bosonic sound waves, is at- 
tractive. Very close to the Feshbach resonance when the 



atomic scattering length is much larger than the effective 
range, the full fcrmion fermion interaction can then be- 
come weakly attractive. Under these circumstances the 
BEC may coexist with a weakly coupled BCS superfluid. 

Over the last two decades many authors, starting with 
Leggett @, H H ES EL studied the BCS-BEC 
crossover in Fermi systems, when the pairing coupling 
constant is varied from very small to relatively large val- 
ues. If the pairing coupling constant is so weak, that 
in vacuum two fermions do not have a bound state, 
one can expect a relatively small pairing gap at a finite 
fermion density, due to the Cooper instability 0] • When 
fcp|a| <C 1 and a < the fermion kinetic energy domi- 
nates over the interaction energy and the Fermi system 
has a positive internal pressure. In the dilute limit (when 
kp\a\ <C 1 and where Uf is the Fermi wave vector and 
a < the scattering length) , the pairing gap is given by 
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Even in the dilute limit the determination of the pre- 
cxponcntial factor requires the evaluation of many-body 
diagrams beyond the naive expected mean-field effects 
alone [l^l IT3 |. At least theoretically, one can contem- 
plate the following question: "How do the properties of 
a Fermi system evolve if one were to increase the pair- 
ing coupling constant until two fermions in vacuum could 
form a bound state?" The arguments presented in Refs. 
00000] lead us to believe that the many fermion 
system will most likely undergo a smooth crossover from 
a weak coupling BCS state to a BEC state, made of rel- 
atively tightly bound fermion pairs. The prospects of 
observing such a crossover were quite dim until the re- 
cent series of experiments 0. While for both posi- 
tive and negative values of the scattering length when 
fcp|a| -C 1 the properties of a many fermion system can 
be described perturbatively, in the intermediate region, 
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when the scattering length is comparable or larger than 
the average inter particle separation the entire system is 
in non-perturbative regime. The basic properties of a 
Fermi system with |<z| = oo have been refined recently in 
a quantum Monte Carlo variational calculation |To| . 

In a Fermi-Bose mixture the interaction responsible 
for the formation of Cooper pairs could be an induced 
interaction 0, 0, 0, 0, 0] . A moving fermion can 
excite a sound wave in the Bose gas, which will perturb 
at a relatively large distance another fermion. This situ- 
ation is similar to the electron-phonon model of BCS su- 
perconductivity. The boson induced interaction between 
fermions is in momentum-frequency representation given 
by 
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where U f b is the fermion-boson coupling constant, q and 
uj are the wave vector and the frequency exchanged be- 
tween the two interacting fermions, e q = h 2 q 2 /2m b , m b 
is the boson mass, n b the boson number density, U bb = 
ATrh 2 a bb /m b and a bb the boson-boson scattering length. 
In the weak coupling limit the two paired fermions ex- 
change a small energy fku — O(sf) and a large linear 
momentum q = 0{kp) and then the induced interaction 
is predominantly static and attractive. It is straightfor- 
ward to show that in coordinate representation the in- 
duced interaction is of Yukawa type [3 [H, E H3, E3 



Ufbf(r) 



II 2 



1 



where 



ft 



U bb An^r 



a bb 



ft 



n b U t 



bb 



2m b s b ^Wnriha^ 



(3) 



(4) 



ft is the healing or coherence length and Sb is the speed 
of the Bogoliubov sound waves in the Bose gas. It is 
notable that the strength of the induced interaction in 
momentum representation at q — is independent of the 
boson density. The radius ft of this interaction could 
exceed significantly abb in a very dilute Bose gas, where 
n b al b <C 1, see Eq. (QJ. 

Using a very simple semiclassical criterion suggested 
by Calogero [19( , one can show that the strength of this 
induced interaction is typically weak. Neglecting for the 
moment the role of £///, Calogero's condition that the 
potential Uf b f{r) can sustain at least one two-particle 
bound state is 
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This criterion is saturated in the case of a simple square- 
well potential and thus cannot be improved. If one as- 
sumes that all scattering lengths and masses are compa- 
rable in magnitude, it is easy to see that this condition 



can hardly be satisfied in the dilute limit, when the heal- 
ing/coherence length is large (ft 3> abb) and the right 
hand side of Eq. © is oc ^Oftb/ft ^ 1- We have thus 
naturally arrived at the conclusion that the net effect of a 
repulsive fermion-fermion interaction £/// and of the at- 
tractive boson induced fermion-fermion interaction Ufbf 
could lead to at most a weak coupling pairing gap in the 
fermion sector. Even in the absence of Uff the induced 
interaction can lead to a weak coupling BCS gap only. 

An extremely interesting system to consider consists 
of fermionic atoms mixed with diatomic molecules of the 
same atom species, as in the recent experiment of Jin's 
group 0. One can change the fermion-fermion scatter- 
ing length to a large and positive value a// := a ^> ro 
(tq is the radius of the interaction, which is also typi- 
cally of the order of the effective range) by means of a 
Feshbach resonance. For the sake of simplicity of the ar- 
gument we shall assume that there are equal amounts of 
two identical fermionic spin subspecies, which we shall 
formally identify with the two spin states of a spin 1/2. 
In such a system one can observe the BCS-BEC crossover 
predicted two decades ago 0, El H , 0, 0. What was 
not appreciated until recently [2fJ is that because not 
only a ff = a > 0, but also a fb = 1-179 a>0pilE^ 
and abb > 0, one can manufacture a metastable mixed 
Fermi-Bose system composed of such fermion atoms and 
weakly bound fermion pairs. The size of the fermion pair 
is 0(a) ^> r and a dilute system of such pairs alone is 
metastable, since a b b = 0.60a > see below and Ref. |2^| 

This is due to the act that the transitions between 
atoms and molecules is a relatively slow process, see Refs. 
0, Q and in particular the discussion of the latest exper- 
iment in Jin's group Q. In fact, at temperatures smaller 
than the binding energy of a pair T <§; h 2 /mfa 2 , the 
pair-pair collision is mostly elastic. Inelastic collisions 
with the formation a tightly bound state and one or two 
fast atoms with kinetic energies of the order ti 2 /m/rg ^> 
h 2 /irifa 2 3> T is possible. The branching ratio for such a 
channel is small (oc (r$/a) a , a~ 3) |2^|. Moreover, since 
the final products will carry a large kinetic energy, they 
will interact with the rest of the slow constituents with 
cross sections of the order r 2 . <C a 2 <C n~ 2/>3 , where n 
is either the fermion or the pair number density. A fast 
atom will not interact with a pair with a cross section of 
the order of a 2 , as one might naively expect, for similar 
reasons why cosmic neutrinos do not have a cross section 
of the order of the Earth radius squared when they im- 
pinge on our planet. Therefore, such process could not 
lead to a noticeable heating of the system. Three-body 
recombination rates can be suppressed by choosing ap- 
propriate low densities. One should recognize however 
that due to the different dependence on the scattering 
length a of the three-body recombination and of the the 
rate of formation of tightly bound molecules, the optimal 
conditions for the formation of an atomic-molecular sys- 
tem correspond to such values of the scattering length 
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for which these rates are approximately equal. The life- 
times of such a mixed atomic-molecular system could be 
estimated as ranging from 10 _2 s to several seconds, de- 
pending on the atom species, peak number density and 
specific value of the scattering length. 

Th total energy density of such a system of spin 1/2- 
fermions is given by the following ex pre ssion (to/ := 
to, TO& := 2m na 3 -C 1, n = n/ + 2nb) 
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where nt = k F /3ir 2 , kp is the Fermi wave vector and 
e-2 = — H 2 /ma 2 is the ground state energy of the pair. The 
specific numerical values quoted here were obtained from 
a detailed calculation of the atom-dimer and dimer-dimer 
scattering lengths described in the next paragraph. There 
are well known corrections of higher order in na 3 -C 1 
to this expression, which we neglect here. The chemi- 
cal potentials for the fermions and bosons in this case 
can be specified independently. Chemical equilibrium 
between fermions and pairs can be established only af- 
ter a very long time, as at these small temperatures 
1^2 1 = fi 2 /ma 2 3> T, see the discussion of recent re- 
sults from Jin's group 4]. The energy density in Eq. 
(|fi*|) has to be supplemented with a contribution arising 
from the weak coupling pairing correlations (condensa- 
tion energy), due to induced interactions, if there is a 
Cooper instability. Since the Fermi subsystem is always 
in the weak coupling BCS limit, this correction to the 
energy is small. 

In the extreme limit when a > ro, where tq is the effec- 
tive range appearing in the low energy expansion of the s— 
wave scattering phase shift k cot 5(k) = —^ + ^r k 2 + ..., 
the fermion-pair scattering length is a/& ~ 1.179 a, a re- 
sult first established by Skorniakov and Ter-Martirosian 
in 1957 Since the value of the atom-atom scatter- 
ing length can be manipulated by means of the Feshbach 
resonance, so can the ratio ro/a, as ro (unlike a) should 
have a small variation with the applied magnetic field in 
the narrow region of the resonance. The dependence of 
a/b as a function of a = a,ff for fixed r is determined, up 
to corrections of order (r /a) 3 , by a simple generalization 
of the Skorniakov and Ter-Martirosian equation |2l|. It 
turns out that the ratio a fb/ a varies by less than 1% 
as ro varies from ro = to ro = 0.3a and the likeli- 
hood of Cooper pair formation is not much affected by 
finite range effects. We have performed as well calcula- 
tions of the dimer-dimer scattering length using a well 
established four-body scattering formalism. The four- 
body equations we solve involve four identical spin 1/2 
fermions (two with spin up and two with spin down), in- 
teracting through a short range potential in 1 Sq alone. 
The formalism we use is identical to the one developed 
in Refs. [25| to study the four-nucleon system at low 



energies and is based on the solution of the Alt, Grass- 
berger and Sandhas (AGS) equations j2(| in the form 
that was proposed by Fonseca and Shanley [13] ■ Since 
this method requires the solution of the underlying two- 
and three-body T-matrix for all relevant channels, it is 
appropriate to mention that we include in the calculation 
three-body sub-amplitudes with total angular momenta 
ranging from l/2 ± to 9/2 ± , since we find that P-wave 
atom-dimer scattering contributes as much as 20 % to 
the dimer-dimer scattering length. The results of the 
calculation indicate that abb = 0.60a. Calculations were 
performed using both a finite and a zero range two-body 
interaction. In the case of a finite range we have varied 
the ratio a/ro from 0(1) to w 2000. The calculations 
of Petrov et al. |23j of the atom-dimer and dimer-dimer 
scattering lengths, based on a novel theoretical approach, 
whose validity we were unable to assess though, agree 
with our result. 

It is remarkable that the interaction properties of this 
mixed system of atoms and weakly bound molecules can 
be described in terms of a single parameter, the fermion 
fermion scattering length a. From the previous analysis 
[lH Hil liH IrH ITq and the arguments given above we 
know now that the tightly bound Cooper pairs, which 
form a BEC, can lead to an additional attraction among 
the fermions. The full fermion-fermion interaction in 
momentum space (and zero frequency) is therefore 
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where hq is the momentum exchanged by two incoming 
fermions with momenta hk\ 2- The s-wave pairing gap is 
determined by the I — partial wave of this interaction 
(if attractive), when both fermion momenta are at the 
Fermi surface (fej — k 2 — kp) 
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otherwise A = 0. The value of the pre-exponential factor 
requires further analysis. The induced interaction has a 
relatively large radius and typically in an extremely di- 
lute system 2fc F £ = 1.127(n / /n 6 ) 1 / 3 (n b a 3 )- 1 / 6 > 1, if 
the densities of the fermion and boson subsystems are 
comparable. Pairing in the s-wave is not likely then, 
since in this case the s-wave full fermion-fermion inter- 
action is repulsive. Still, for appropriate values of the 
parameters s-wave pairing is possible, see Fig. ^ F° r 
other values of the parameters p-wave pairing should be 
considered |28j , but the size of the p-wave gap will typi- 
cally be smaller. 

The spectrum of sound waves in such a system has 
two branches. The bosonic pairs will naturally have a 
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FIG. 1: The ratio A/sf (sf = h 2 k F /2m) as a function of 
rhf/rhb, for a fixed boson number density rib = 10 13 cm -3 
and riba 3 = 0.064 (solid line) and n^a 3 = 0.037 (dashed line) 
respectively. The dots show the value of the gap in the case 
of p-wave paring for rn,a i = 0.064. 

Bogoliubov sound branch with the speed Sh, see Eq. 
while in the fermion sector the sound waves will have 
the speed :29j, s f = v F /V3 = h^rif) 1 / 3 /%/3m, un- 
like in a normal Fermi gas, where the speed of the Lan- 
dau's zero sound waves is larger then vf 30]. The ra- 
tio of these two sound speeds could in principle have 
any value, depending on the specific values of the scat- 
tering length and fermion and boson number densities 
Sf/s b oc (ny/n fc ) 1 / 3 /(n fc a 3 ) 1 / 6 . For dilute systems with 
comparable amounts of bosons and fermions the slowest 
mode is the bosonic Bogoliubov sound wave. There are 
relatively small corrections to the values of Sb and s/, 
due to interactions and also due to some mixing between 
the two types of sound waves. We simply mention that 
for these systems (in both normal and superfluid fermion 
phases) there is no long-wave instability of the type sug- 
gested in Ref. 31]. Even though the present analysis 
was limited to the homogeneous case, the properties of 
a system in a trap can be easily inferred and as, it was 
shown in Ref. |3^ | , an atom- molecule mixture of the kind 
considered here acquires a complex and quite unexpected 
spatial structure. 

Discussions with G.F. Bertsch, D.S. Petrov, G.V. 
Shlyapnikov, B. Spivak, S. Stringari and L. Viverit 
are greatly appreciated. AB and PFB acknowledge 
DOE support and ACF was supported in part by grant 
POCTI/FNU/37280/2001. 
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